Resolving the Host Galaxy of the Nearby QSO I Zw 1 with 
Sub-Arcsecond Mult i- Transit ion Molecular Line Observations 
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ABSTRACT 

We present the first sub-kpc (~ 0.7" s» 850 pc) resolution 12 CO(1-0) molecular line obser- 
vations of the ISM in the host galaxy of the QSO I Zw 1. The observations were obtained 
with the BIMA mm-interferometer in its compact A configuration. The BIMA data are comple- 
mented by new observations of the 12 CO(2-I) and 13 CO(I-0) line with IRAM Plateau de Bure 
mm-interferometer (PdBI) at 0.9" and 1.9" resolution, respectively. These measurements, which 
are part of a multi- wavelength study of the host galaxy of I Zw 1, are aimed at comparing the 
ISM properties of a QSO host with those of nearby galaxies as well as to obtain constraints on 
galaxy formation/evolution models. Our images of the 12 CO(T0) line emission show a ring-like 
structure in the circumnuclear molecular gas distribution with an inner radius of about 1.2 kpc. 
The presence of such a molecular gas ring was predicted from earlier lower angular resolution 
PdBI 12 CO(1-0) observations. A comparison of the BIMA data with IRAM PdBI 12 CO(2-l) 
observations shows variations in the excitation conditions of the molecular gas in the innermost 
1.5" comprising the nuclear region of I Zw 1. The observed properties of the molecular cloud 
complexes in the disk of the host galaxy suggest that they can be the sites of massive circumnu- 
clear star formation, and show no indications of excitation by the nuclear AGN. This all indicates 
that the molecular gas in a QSO host galaxy has similar properties to the gas observed in nearby 
low luminosity AGNs. 

Subject headings: galaxies: ISM, nuclei — quasars: individual (I Zw 1) 
1. Introduction 
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One of the open questions in galaxy evolution 
is concerned with the formation of active galac- 
tic nuclei (AGN) and its relation to star forma- 
tion (SF) in early type galaxies. The separation 
of starburst and AGN components in extragalac- 
tic objects - especially in host galaxies of quasars 
and QSOs - is a key problem in the investiga- 
tion of evolutionary sequences proposed for AGNs 
(Norman & Scoville 1988; Sanders et al. 1988; 
Rieke et al. 1988; Haas et al. 2003). Although it 
is not exactly known how the host galaxy affects 
the energy release of the QSO, there is statistical 
evidence for preceding mergers or current interac- 
tion of QSO host galaxies with a companion galaxy 
(McLeod k Rieke 1994;Lim&Ho 1999). It is be- 
lieved that quasar activity is a common, but short 
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lived, phenomenon in galaxy evolution (McLeod 
ct al. 1999). Furthermore, host galaxies of QSOs 
show enhanced SF activities (Courvoisier 1998). 
The transition between Ultra Luminous Infrared 
Galaxies (ULIRG) and AGN seems to be contin- 
uous, all showing signs of enhanced star forma- 
tion in their nuclear regions (Genzel et al. 1998; 
McLeod et al. 1999). 

Millimeter molecular line observations are ideal 
for studying the mechanisms which transport 
the molecular gas into the AGN. Observations 
of abundances, excitation and dynamics of the 
molecular interstellar medium in the central re- 
gions of AGN are essential, since the interstellar 
matter provides the "fuel" for star formation as 
well as the central engine. Consequently, there 
are large numbers of high resolution observations 
available for nearby objects. These observations 
have revealed the presence of circumnuclear star- 
burst rings in a large number of (nearby) active 
and IR luminous galaxies (e.g. NGC 1068, Plane- 
sas et al. 1991; NGC7469, Genzel et al. 1995; 
BIMA SoNG, Heifer et al. 2003). The extension 
of high angular resolution observations of molecu- 
lar gas emission lines to QSO hosts is imperative 
in order to understand the connection between 
local active galaxies and high-z QSOs. The re- 
sults can be used to refine model predictions of 
the physical conditions in high z- QSOs (such as 
used e.g. in Combes et al. 1999). 

Due to the limited angular resolution of sin- 
gle dish mm-wavelength telescopes, only interfer- 
omctric observations at these wavelengths allow 
insight into the morphology and into the kinemat- 
ics of molecular clouds in the nuclear region of 
QSO host galaxies with sufficient angular resolu- 
tion. High resolution molecular line observations 
allow for detailed kinematic studies of the cold in- 
terstellar medium and the derivation of important 
parameters such as gas masses, surface mass den- 
sities and star formation efficiencies. Observations 
of multiple CO transitions reveal excitation condi- 
tions and thus can be used to constrain the physi- 
cal conditions in an observed region. They can be 
used to better constrain the contribution of star 
formation to the total observed infrared luminosi- 
ties. 



1.1. The nearby QSO I Zw 1 

The radio-quiet QSO I Zw 1 is regarded as the 
closest QSO which can be used for detailed studies 
of its host galaxy (Tab. 1). I Zw 1 has a systemic 
velocity of 18,290 km s , which corresponds to a 
redshift of 0.0611 (Condon et al. 1985), or a dis- 
tance of 255 Mpc (h = 0.72, Spergel et al. 2003). 
The nucleus of I Zw 1 is extremely bright in the op- 
tical (M B of -23.45 mag, Schmidt & Green 1983; 
Barvainis et al. 1989)) and also has very bright 
X-ray emission (Kruper et al. 1990). I Zw 1 
belongs to the class of infrared luminous galax- 
ies (L FIR = lO 11 - 3 ^©, Haas et al. 2003), how- 
ever its QSO nature is not apparent in the FIR. 
There is a strong indication for interaction be- 
tween I Zw 1 and a neighboring companion (e.g. 
Lim k Ho 1999; Scharwachter ct al. 2003). 

A direct comparison between the FIR emission 
and CO molecular line emission in the host galaxy 
reveals that the FIR continuum emission of I Zw 1 
is predominantly thermal in nature (Barvainis et 
al. 1989). The FIR to 12 CO(1-0) ratio shows a 
star formation efficiency of about 30 Lq/Mq in 
the nuclear region (Eckart et al. 1994). Schin- 
ncrcr ct al. (1998) show that the molecular gas 
mass in the central 3.3" (3.7 kpc) is 2/3 of the total 
observed mass of 9 x 1O 9 M . A kinematic analy- 
sis of their 1.9" resolution PdBI 12 CO (1-0) obser- 
vations suggests that the central molecular gas is 
distributed in a ring of 1.5" (~ 1.9 kpc) diameter. 
However the observations do not spatially resolve 
this ring. High resolution NIR observations show 
that the nuclear NIR spectrum is dominated by 
emission from the AGN, nevertheless about 25% 
of the NIR emission can be attributed to a nu- 
clear starburst (Schinnerer et al. 1998), which 
might be located in the molecular gas ring. The 
star formation activity in the host galaxy of I Zw 1 
underlines the important role of starbursts in the 
evolution of QSOs. 

Only sub-arcsecond molecular line observations 
would allow actual imaging of this possible star- 
burst ring. The extension of BIMA's longest base- 
lines up to a maximum of almost 2 km makes the 
array the millimeter instrument with the highest 
angular resolution currently available. Here we 
present 12 CO(1-0) ~ 0.7" angular resolution ob- 
servations of the QSO I Zw 1. This corresponds 
to a spatial resolution of roughly 850 pc. We also 
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present new observations with the IRAM PdBI in 
12 CO(2-l) and 13 CO(1-0) line at 0".9 and 2."1 res- 
olution, respectively. 

2. Observations and data reduction 
2.1. BIMA observations 

The data presented here were taken in Decem- 
ber 1999 in the compact A array configuration 
(A-) of the ten element BIMA array (Welch et 
al. 1996) with baselines between 80 and 1310 m. 
Four more tracks of A array observations were ob- 
tained in February 2000, December 2001, and Jan- 
uary 2002, but all of those yielded only poor qual- 
ity data, demonstrating that this project poses 
requirements which are at the limit of what can 
be done with long baseline mm-interferomctry to- 
date. The digital correlator was configured to 
cover two 800 MHz bands with the upper side- 
band centered at the redshiftcd 12 CO(1-0) fre- 
quency of 108.633 GHz. The sky opacity is signifi- 
cantly lower at the frequency of redshifted 12 CO(l- 
0) transition than at the rest frequency. System 
temperatures during the observations ranged from 
180 to 400 K, single-sideband. We flagged data 
such that most observations with system temper- 
atures above 300 K were rejected. The data were 
reduced with the MIRIAD software package (Sault 
et al. 1995). 

The fluctuations in the atmosphere follow a 
Kolmogorov power law distribution up to the 
maximum baseline of 2 km in the A array con- 
figuration of the BIMA array. This means that 
phase variations increase with baseline length, and 
therefore long baseline observations require a par- 
ticularly careful phase calibration. We employed 
fast phase referencing: the observations were 
switched between the source, the phase calibra- 
tor (0108+015, distance from I Zw 1 - 11°) and 
an additional point source, the quasar 0121+118 
(angular distance from I Zw 1 ~ 7°) which serves 
as a test source to determine the accuracy of the 
solutions for the phases. One complete cycle takes 
less than two minutes in order to follow closely the 
atmospheric phase (Holdaway & Owens 1993). 

The test calibrator was used to select the best 
quality data: only observing cycles for which the 
test source maps to a point source down to the 
20% level of its peak flux were selected (Fig. 1). 
This criterion was chosen to allow for a good qual- 



ity image of the source. The reverse check, us- 
ing 0121+118 as phase calibrator and 0108+015 
as the test source, verified that the phase solu- 
tions derived from the test calibrator are of the 
same accuracy. Only the data taken in Decem- 
ber 1999 have sufficient phase coherence to yield 
sub-arcsecond resolution observations. The data 
presented here are only from this observing track. 
The phase solutions derived from the test source 
0121+118 were used for the images shown here, 
since this source is closer to I Zw 1 and therefore 
provides the best phase coherence. 0121+118 has 
a very flat spectrum with a flux of 1.1 Jy at 3.7 cm, 
2 cm, and 7 mm (VLA calibrator list) and we use 
the same value for our 3 mm band. We estimate 
the uncertainty in the amplitude calibration to be 
better than 30%. In order to increase the signal 
to noise ratio in the data, we smoothed them to 
an effective angular resolution of 0.7". 

2.2. IRAM PdBI observations 

I Zw 1 was observed simultaneously in the 
12 CO(2-l) and 13 CO(1-0) line in August 1996 with 
the four-clement IRAM Plateau de Bure millime- 
ter interferometer (PdBI) in the C configuration. 
Further dual-frequency observations in the B con- 
figuration were obtained in January 1997 and Jan- 
uary/February 1998 with the then five-element 
array. The correlator was centered on the red- 
shifted 12 CO(2-l) frequency of 217.26321GHz at 
1 mm and the redshifted 13 CO(1-0) frequency of 
103.85578 GHz at 3 mm. With baselines between 
40 and 280 m the data have a spatial resolution of 
~ 0.86" with uniform weighting for the 12 CO(2-l) 
line (Fig. 2), and ~ 2.1" with natural weighting 
for the 13 CO(1-0) line. The noise per 20(40) km 
s _1 wide channel is 1.7(0.33) mJy beam -1 in the 
combined data for the 12 CO(2-1)( 13 CO(1-0)) line. 
The quasar 3C454.3 served as a passband calibra- 
tor, while standard IRAM PdBI phase calibrators 
were observed every 20 minutes. The data were 
calibrated and mapped using the IRAM GILDAS 
software package (Guilloteau & Lucas 2000). 

3. Results 

3.1. Spatial distribution of the 12 CO(1-0) 
line emission 

The gas disk of I Zw 1 is clearly resolved in its 
12 CO(1-0) line emission at 0.7" resolution (Fig. 
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Table 1 

Adopted Properties of I Zw 1 (=PG0050+124) 



Property 


Value 


Right ascension (J2000.0) 


00 h 53 m 34 s .9 


Declination (J2000.0) 


+12°41'36".2 


Inclination (cleg) 


38 


Position angle (deg) 


135 


Systemic velocity (km s _1 ) 


18,290 


Distance (Mpc) 


255 


1" 


1.24 kpc 



~The coordinates refer to the phase center 
of the millimeter observations. The inclination 
and position angle are taken from Schinnerer 
ct al. (1998), while the systemic velocity is 
from Condon et al. (1985). 



Table 2 

i zw 1 12 co(1-0) beam averaged parameters 



Region 


Av 


T B 


^co 


7V H2 




[km s _1 ] 


[K] 


[K km s- 1 ] 


[10 21 cm" 2 ] 


Disk (9.0") a 


310 


0.03 


10 


2 


Nucleus (2.0") 


400 


1.8 


720 


144 


kpc ring (0.7") 


100 


12 


1200 


240 



a from Schinnerer et al. (1998) 



4 



3) . These are the first molecular line observations 
which resolve the emission in the nuclear region 
and show that the emission is not centrally peaked: 
the strongest molecular line emission is not coin- 
cident with the optical center of the galaxy, but 
situated at a distance between 1" and 3", mainly 
to the South and East of the galaxy center. This 
corresponds to a distance of between 1.2 kpc and 
3.7 kpc from the nucleus. This emission, which 
occurs at positive offsets in R.A. with respect to 
the major axis, has predominantly positive veloc- 
ities (see position- velocity (pv) diagram in Figure 

4) . When the velocity integrated BIMA 12 CO(l- 
0) data are smoothed to a resolution of 2" (Fig. 5, 
top left), which corresponds to the resolution of 
the corresponding 12 CO (1-0) PdBI observations 
(Fig. 5, top right), we find a positional offset of 
~ 1" in the center of emission of both images with 
respect to each other. The measured peak flux in 
the BIMA image, which is smoothed to a resolu- 
tion of 2" (Fig. 5, top left), is 28 mJy beam -1 
or 11 Jy beam -1 km - s, slightly larger than the 
12 CO(1-0) peak flux of 8 Jy beam" 1 km -1 s, at the 
same angular resolution observed with the PdBI. 
The difference of ~ 38% is still within the expected 
range, if we assume an uncertainty of the PdBI 
calibration of 20%. 

This significant spatial offset of the 12 CO(1-0) 
center of emission in the BIMA and PdBI obser- 
vations can not be due to a pointing error of the 
BIMA observations, since a pointing error of this 
magnitude is not visible in the image of the test 
quasar (Fig. 1). The fact that the negative velocity 
components in the BIMA observations (Fig. 4) are 
less pronounced than in the lower resolution PdBI 
data (Fig. 6, left) indicates that the molecular 
emission at these velocities has a more extended 
distribution and therefore is partially resolved out 
in the high resolution BIMA data. This effect nat- 
urally explains the shift of the observed center in 
the BIMA data. Probably the most convincing 
indication that the denser, and therefore probably 
more spatially confined, molecular clouds are situ- 
ated to the South-East of the nucleus comes from 
the 13 CO(1-0) observations we obtained with the 
PdBI (Fig. 5, bottom left): the center of emission 
in these data is also shifted towards the South- 
East. It is worthwile to note that the peak of 
the line emission in all CO transitions is shifted 
by ~ l"to the South-East of the nominal optical 



center. This could be due to the fact that the dis- 
tribution of the CO line intensity is asymmetric 
with respect to the AGN. This assumption could 
be tested if the position of the dynamical center 
could be accurately determined; however, even the 
resolution of our 12 CO(2-l) data is not adequate 
to do so. 

3.2. Properties of the molecular gas 

The observed pv diagram along the major axis 
at 0.7" resolution (Fig. 4) shows a circular ro- 
tation which is consistent with the lower resolu- 
tion PdBI 12 CO(1-0) observations (Fig. 6, left; 
see also Scharwachter et al. 2003). The con- 
sistency of the observed p — v diagrams in the 
BIMA and PdBI observations strongly supports 
the interpretation that the spatial offset of the 
peak emission in the BIMA image is not an ar- 
tifact due to calibration errors. The peak emis- 
sion of 35 mJy beam -1 in the BIMA pv diagram, 
which is observed over a velocity interval of about 
100 km s -1 , corresponds to an average bright- 
ness temperature of Tg = 12 K in this velocity 
interval. This indicates a large beam filling fac- 
tor, respectively a significant amount of molecular 
gas within a region of 600 pc 2 . Schinnerer et al. 
(1998) show that the standard conversion factor 
iV(H 2 )//(CO) = 2xl0 20 cm -2 K -1 km -1 scan be 
applied to the observed 12 CO (1-0) line emission 
from I Zw 1. 

The peak flux of 18 mJy beam -1 in the ve- 
locity integrated 0.7" resolution image (Fig. 3) 
translates into a molecular hydrogen column den- 
sity of 4.9 x 10 23 cm -2 . About 50% of this emission 
has a velocity consistent with its origin being in a 
circumnuclear ring at a distance of ~ 1.5". The 
corresponding column density of 2.4 x 10 23 cm~ 2 
in the molecular ring is two orders of magnitude 
larger than the average molecular gas column den- 
sity of 2 x 10 21 cm~ 2 in the disk region, and roughly 
a factor of two larger than the average column den- 
sity of 1.4 x 10 23 cm~ 2 in the nuclear region (the 
innermost 2" region: see Table 2 for a summary of 
the observed 12 CO(1-0) molecular line properties). 

3.3. Multi-transition analysis of the molec- 
ular line emission 

Our PdBI observations of the 12 CO(2-l) and 
13 CO(1-0) line detect an integrated line flux lco 
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of 18Jykm s 1 and 0.5Jykm s , respectively. 
Thus they recover only about 30% (15-20%) of the 
observed IRAM 30m single dish flux of 61Jykm 
s- 1 (3.1Jykm s" 1 ) for the 12 CO(2-l) ( 13 CO(l- 
0)) line (Eckart et al. 1994). In contrast to the 
BIMA 12 CO (1-0) observations, the 0.86" resolu- 
tion PdBI 12 CO(2-l) observations (Fig. 5, bottom 
right) do not show any particular spatial structure 
in the nuclear region. This discrepancy can only 
be explained by a difference in the excitation of 
the molecular gas in the nucleus (innermost 1") 
and the emission from molecular clouds which are 
situated further out. This interpretation is sup- 
ported by both the PdBI and BIMA 12 CO(1-0) 
observations: A comparison of the pv diagrams 
from the PdBI 12 CO(1-0) observations from Schin- 
nerer et al. (1998) (here chosen for the better 
signal to noise) and our new PdBI 12 CO(2-l) and 
13 CO(1-0) observations can be seen in Figure 6. 
The peaks of the line intensity are shifted from 
±150 kms -1 in the 12 CO(1-0) line to lower veloc- 
ities of ±(50 - 100) kms" 1 in the 12 CO(2-l) line. 
This difference becomes even more obvious in the 
^ >21 ~ i2co(i-o) nne P e& k ratio map of the major 
kinematic axis pv diagram shown in Figure 7. For 
these ratios the resolution of the 12 CO(2-l) obser- 
vations was smoothed to match the PdBI 12 CO(l- 
0) observations. The line peak brightness ratio in 
the inner <~ 1.5" of I Zw 1 is « 0.5. This drops 
by almost a factor of two to w 0.3 at the position 
of the 12 CO(1-0) line peaks at ±150 km s" 1 . This 
suggests that the outer 12 CO(1-0) clouds, which 
are apparently not resolved out by the BIMA ob- 
servations with a similar angular resolution as the 
PdBI 12 CO(2-l) observations, are more subther- 
mally excited than the molecular gas in the inner- 
most region of I Zw 1. The fact that the denser 
and more thermalized molecular material breaks 
up into several distinct velocity components in the 
12 CO(2-l) emission (Figure 6) strongly suggests 
that circumnuclear star formation and not the cen- 
tral AGN is exciting the molecular line emission. 

The derived r w = iii co i_ line peak bright- 
ness ratio for the three 13 CO(1-0) line peaks above 
3a (see Figure 6, right) is about 15. This value can 
probably be regarded only as an upper limit due to 
the large amount of missing flux which has been 
resolved out by our interferometric observations, 
in particular in the 13 CO (1-0) observations (see 
above). We therefore avoid further interpretation 



of these ratios. 

NIR spectra of the inner 3.3" of I Zw 1 indi- 
cate the presence of circumnuclear star formation 
(Schinnerer et al. 1998). The millimeter line ra- 
tios for the inner 1.5" suggest that these GMCs 
could be the site of this massive star formation. 
However, the BIMA observed 12 CO(1-0) GMCs, 
which are further out, are also likely candidates 
for the observed starburst due to their high col- 
umn densities. A starburst in these clouds might 
be triggered by the interaction of I Zw 1 with the 
companion galaxy (see Scharwachter et al. 2003). 

3.4. Molecular gas kinematics 

There is some enhanced emission in the BIMA 
pv diagram (Fig. 4) at velocities which are "forbid- 
den" with respect to circular motion, most promi- 
nently at a positive angular offset of 2" along the 
major axis between —40 and —60 kms -1 . Other 
than the emission from the circular rotation, it is 
pronounced only in one bin of the p-v diagram, 
and with a signal to noise ratio of 4cr its signif- 
icance is not sufficient to be considered a clear 
detection. The 12 CO(2-l) velocity field (Fig. 8) 
also shows indications of deviation from pure cir- 
cular rotation for radii of r > 1", as is evident in 
the change in the line of nodes. Due to the lack 
of emission, this potential non-circular motion can 
not be followed to larger radii. Together with the 
"forbidden" 12 CO(1-0) it can be seen as circum- 
stantial evidence that some transportation mech- 
anism for the molecular gas is active in the disk of 
I Zw 1. However, data of higher angular resolution 
and improved sensitivity will be needed to confirm 
this. A more detailed analysis of the observed ro- 
tational pattern can be found in Scharwachter et 
al. (2003). 

4. Summary and Conclusions 

Our BIMA 12 CO(1-0) observations represent 
the first molecular line observations of the host 
galaxy of a QSO with sub-kpc resolution. The 
molecular gas disk of a QSO is for the first time 
spatially resolved into individual giant molecu- 
lar cloud complexes. We observe giant molecular 
clouds with peak molecular hydrogen column den- 
sities of 2.4 x 10 23 cm~ 2 at a distance between 1" 
and 3" from the nucleus. With such high column 
densities these clouds could be actively forming 
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stars. 

The combination of BIMA observations with 
new and previously published PdBI CO data al- 
lows for the first interfcrometric multi-transition 
study in the disk of I Zw 1. The line brightness 
ratios in the inner 1.5" are consistent with moder- 
ately dense cold GMCs, and they are not peaked at 
the center. This strongly suggests that the AGN 
has no significant effect on the central molecular 
material. 

The exact location of the circumnuclear star- 
burst very likely seen in nuclear NIR spectra 
(Schinnerer et al. 1998) is still not identified. Pos- 
sible candidates are the resolved 12 CO (1-0) molec- 
ular clouds seen by BIMA as well as the more ther- 
malized material traced by the higher line ratios 
inside the inner 1.5". However, any of these molec- 
ular clouds are likely sites of a massive starburst 
which contributes to the observed far-infrared lu- 
minosity of L F i R — 10 113 LQ. 

The distribution and properties of the molec- 
ular gas in the host galaxy of the nearby QSO 
I Zw 1 are quite similar to what is observed in 
nearby low luminosity AGNs (e.g. Paglionc et al. 
2001; NUGA Survey, Garca-Burillo et al. 2003). 
This suggests that the ISM in QSO at high rcdshift 
might be similar to nearby low luminosity AGN as 
well. 

Four additional attempts to re-observe I Zw 1 
in 12 CO (1-0) with BIMA in the A configuration 
all failed to provide useful data, due to insufficient 
phase coherence of the observations. This demon- 
strates that the observations presented here are 
on the limit of what can be done without active 
phase correction systems, at least from a site such 
as Hat Creek where the BIMA array is situated. 
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constructive comments from the anonymous ref- 
eree are much appreciated. A. Eckart. and J. 
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494. J. Scharwachter is also supported by a schol- 
arship for doctoral students of the Studienstiftung 
des deutschen Volkes. 
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Fig. 1— BIMA image of the test source 0121+118 
in 12 CO(1-0). The contour lines start at 20% of 
the peak level in increments of 10%. Only the 20% 
contour shows a slight deviation from the ideal 
point source response. The synthesized beam is 
shown in the lower left of the image. 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.0. 
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Fig. 2. — Channel maps of 12 CO(2-l) line emission 
in I Zw 1 observed with the PdBI. The relative 
velocities with respect to the observed systemic 
velocity are indicated in the upper left corner of 
each panel. The data have a resolution of 0.86" 
and the contours are in steps of 5mJy beam -1 
(~ 3<r). The rotation of the galactic CO disk is 
clearly evident in the channel maps. The clean 
beam is shown for reference in the lower left corner 
of each panel. 
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Fig. 3.— BIMA image of the 12 CO(1-0) emis- 
sion smoothed to a resolution of 0.7". The 
solid contours show the velocity integrated emis- 
sion between —160 and 240 km s -1 in la steps 
(3mJy beam -1 ) starting at 4 a (9 mJy beam -1 ). 
If present, dashed contours would represent the 
corresponding negative values for the flux density. 
The ellipse in the center of the image marks the 
projected distance of 1 kpc (0.8") from the center. 
The diagonal line indicates the major kinematic 
axis of the galaxy. The ellipse in the lower left 
corner shows the synthesized beam. 
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Fig. 4.— BIMA 12 CO(1-0) position-velocity dia- 
gram along the major axis (135°) of I Zw 1 with 
a spatial resolution of 0.7". The contours start at 
60% (3 a) of the peak flux of 35 mJy beam -1 with 
increments of 15% (lcr). 
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Fig. 5.— Top left: BIMA image of the 12 CO(1-0) 
emission smoothed to a resolution of 2.0". The 
solid contours show the velocity integrated emis- 
sion between —160 and 240 km s -1 in lcr steps 
(4 mJy beam -1 ) starting at 4 a (17mJy beam -1 ). 
Dashed contours would show the corresponding 
negative values for the flux density, if it were 
present in the image. The ellipse in the lower 
left corner shows the synthesized beam. To ob- 
tain integrated line fluxes, these and the follow- 
ing values should be multiplied by the velocity 
range. Top right: PdBI 12 CO (1-0) Oth moment 
image with 1.9" resolution, integrated between 
-240 and +240 km s -1 . The contours start at 10% 
of the peak flux of 0.016 Jy beam -1 and arc in 
steps of 10%. Bottom left: PdBI 13 CO(1-0) im- 
age integrated over the velocity range of —240 and 
+240 km s -1 . The contours are at 0.3 (3a), 0.4, 
0.5, 0.7, and 0.8 mJy beam -1 (7a). Bottom right: 
PdBI 12 CO(2-l) Oth moment image at 0.86" reso- 
lution, integrated between —240 and +240 kms -1 . 
The contours start at 10% of the peak flux of 0.017 
Jy beam -1 and are in steps of 10%. A cross marks 
the position of the optical center in each panel. 
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Fig. 6. — PdBI data pv diagrams along the ma- 
jor axis (135°) of I Zw 1. The 12 CO(1-0) emis- 
sion (left) is shown at 1.9" resolution, the 12 CO(2- 
1) emission (middle) at 0.86" resolution and the 
13 CO(1-0) line emission (right) at 2.1" resolution. 
The contours start at 3a in steps of 3<7 for the 
12 CO data with la = 2(1.7) mJy beam" 1 for 
the 12 CO(1-0) ( 12 CO(2-l)) line. For the 13 CO(l- 
0) the contours start at 2a in steps of la with 
la = 0.3 mJy beam -1 . 
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Fig. 8.— PdBI 12 CO(2-l) velocity field at 0.86" 
resolution. The velocity field is shown in contours 
(the thick line represents Vhel — 18,316 km s -1 , 
the solid (dashed) lines are positive (negative) ve- 
locities in steps of 20km s _1 ). The 0th moment 
image is shown in grayscale for comparison. The 
beam is indicated in the lower left corner. 



Fig. 7.— PdBI 12 CO line pv diagrams at 2.1" 
resolution. The i?2i = i^co(i-o) P ea ^ brightness 
ratio (middle) clearly shows a higher ratio in the 
inner 1.5" which is decreasing by almost a factor 
of two towards higher velocities of ±150 km s -1 . 
The contours are at 0.3 (black line), 0.4 (grey line) 
and 0.5 (white line). For reference the 12 CO(1-0) 
(left) and 12 CO (2-1) (right) pv diagrams are shown 
at a resolution of 2.1". The contours start at 3a 
in steps of 3a with la of 1.8 mJy beam -1 . 
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